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Ethanol injection and variations of it are a class of methods where two miscible phases — one of which contains dissolved lipids — 
are mixed together leading to the self-assembly of lipid molecules to form liposomes. This method has been suggested, among 
other applications, for in-situ synthesis of liposomes as drug delivery capsules. However, the mechanism that leads to a specific 
size selection of the liposomes in solution based self-assembly in general, and in flow-focussing microfluidic devices in particular, 
has so far not been established. Here we report two aspects of this problem. A simple and easily fabricated device for synthesis of 
monodisperse unilamellar liposomes in a co-axial flow-focussing microfluidic geometry is presented. We also show that the size 
of liposomes is dependent on the extent of micro-convective mixing of the two miscible phases. Here, a viscosity stratification 
induced hydrodynamic instability leads to a gentle micro-mixing which results in larger liposome size than when the streams are 
mixed turbulently. The results are in sharp contrast to a purely diffusive mixing in macroscopic laminar flow that was believed to 
occur under these conditions. Further precise quantification of the mixing characteristics should provide the insights to develop 
a general theory for size selection for the class of ethanol injection methods. This will also lay grounds for obtaining empirical 
evidence that will enable better control of liposome sizes and for designing drug encapsulation and delivery devices. 
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1. Introduction 

Liposomes offer a great potential in drug delivery systems 
( |Lasic|[T992] ) in terms of versatility in size and charge, ability to 
encapsulate hydrophilic (inside the aqueous core) and lipophilic 
drugs (within the lipid bilayer), relative non-toxicity compared 
to other carrier systems, ability to escape macrophages in the 
lungs (stealth liposomes) and ease of preparing covalently cou- 



2002| ) . Ethanol and water being miscible in all proportions, 



leads to a quick mixing of the phases resulting in induction of 
the lipid self-assembly into bilayers. The initial structures are 
disc-like membranes, which grow and close upon themselves 
or fuse with other discs to form vesicular compartments. The 
size of these lipid vesicles ranges from about 20 nm up to about 
a few micrometers and they may be composed of one or more 
concentric membranes (unilamellar and multilamellar vesicles), 
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2005] |Lasic| |1998| ). Liposome size and size distribution are 
critical parameters to control, since these influence drug dosage 
and the clearance rate of the drug from the body. Liposome 
sizes in the range 100 nm to 200 nm is most suitable for passive 
targeting of cancerous tissue by enhanced permission and reten- 
tion mechanism (EPR) (Mae da et al.||2000t|Lembo and Cavalli| 
|20T0l ). 

One of the ways liposome formation occurs is when of lipids 
molecules self-assemble from a solution phase due to changed 
solvent conditions, such as increase in hydrophobic costs due 
to the entry of water. This phenomenon may be classified un- 



der a broad principle of coacervation pshii et aLj |1995| ) and 
ethanol injection is one such coacervation methods in which 
lipids dissolved in ethanol are brought into contact with wa- 
ter ( [Batzri and Korn 1973[ Kremer et al. 1977 Wagner et al 



There is only a general understanding of the phenomenon of 
liposome formation in spontaneous self-assembly — formation 
of discs of the bilayer membrane ( Lasic| |1992| ) which close 
upon to form vesicles after a critical size when the edge energy 
balances the bending energy (Helfrich| 1973| ). There has so 
far been no definite prediction of the size of liposomes formed 
in the ethanol injection methods, mainly owing to method of 
mixing the two phases and a wide distribution of sizes. The 
size prediction and control, as outlined above, is crucial from a 
drug-delivery point of view. 

In the general class of ethanol injection methods suggested 
in the literature is the microfluidic method, offering a fine con- 
trol on manipulation of fluid and interfaces, and in small vol- 



umes of fluids in the picoliter range. Recently, Jahn et al. ( 2004 
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2007 ) introduced a microfluidic two stream hydrodynamic-focussin^ 
method ( [Qttino and Wiggins, 2004) which was shown to gener- 
ate liposomes of varying sizes and distributions by controlling 
the flow rates. In this device, phospholipids dissolved in a water 
miscible solvent (such as ethanol) is sent through central mi- 
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crochannel and water through a side microchannel in a planar 
geometry. 



Jahnet al. (2007) suggested that an interplay between molec- 



ular diffusion of alcohol and laminar convection results in the 
formation of liposomes of controllable distribution. The size 
control was achieved in the range 50 to 150 nm by adjusting the 
ratio of flow rates and the over all flow rate. While the obser- 
vations provided an empirical way to achieve a control of the 
liposome size, we do not have a quantitative model for predict- 
ing the sizes. This lack of knowledge does not confine itself to 
the micro fluidic method alone. In the entire gamut of ethanol 
injection methods ([Batzri and Korn||1973[|Kremer et al.|[T977} 



Wagner et al.| |2002[ |Ishii et aLj |1995| ) there is no theoretical 



model yet that can explain the size selection or the formation of 
multi-lamellar liposomes in some methods. 

We hypothesise that the nature of flow and mixing has a sig- 
nificant effect on determining the outcome of the self-assembly. 
To show this, we first choose a simpler configuration — core- 
annular two stream flow, and compare the liposome formation 
in this system along with other methods of mixing the two 
streams. We find interesting connections of the liposome size 
with the nature of micro-mixing. 

An additional contribution in this paper is in the method of 
device fabrication itself. While the channel geometry suggested 
in |Jahn et aL] ( |2004| ) is easy to replicate, it requires sophisti- 
cated instruments to create the initial moulds. The channel flow 
also poses a theoretical difficulty in understanding the interac- 
tion of flow and diffusion, as there are no simple analytical so- 
lutions to the velocity and concentration profiles. A co-axial 
core-annular (or core-sheath) flow is easily amenable to simple 
calculations owing to the axial symmetry of the solution. Here 
we introduce a simple method to fabricate a co-axial hydrody- 
namic focussing device to synthesise liposomes following the 
ideas of glass-PDMS hybrid microfluidic device ( Jeong et al 



2004), and a method of swelling of PDMS in an organic sol- 



vent to create channels (Ver ma et al.||2006| ). All of these can 
be easily reproduced in any laboratory without requirement of 
sophisticated instruments and environments to carry out photo- 
lithography (such as, spin-coating, mask alignment tools, gold 
room etc.). 

We first outline the procedure to fabricate the co-axial (core- 
annular) flow generating device in the Methods section. Fol- 
lowing this, we show various studies using the axial flow-focussing 
method — the effect of flow rates, type of mixing, lipid concen- 
tration, ethanol concentration, lipid chemistry on the size of li- 
posomes. We then present evidence that supports our claim that 
it is micromixing that leads to selection of larger size liposomes 
formed within the micro channel. 

2. Materials and Methods 

2.1. Materials 



(18: 1-DOPC) (Sigma Aldrich), and Hydrogenated Soy-phosphocholine 
(gift sample from Lipoid) is used as purchased without further 
purification. Ethanol (AR grade, Merck) is used to dissolve 
phospholipids. 

All the lipids from Avanti Polar Lipids were in a solution 
form with chloroform (20 mg per ml) as solvent. An appro- 
priate volume of lipid solution is taken out in a round bottom 
flask (RBF) and chloroform is evaporated from lipid solution 
using a rotary evaporator. This forms a thin film on the glass 
surface of the RBF. The RBF is then kept in a vacuum desic- 
cator overnight to remove traces of chloroform. An appropriate 
volume of ethanol is added to the RBF to make various concen- 
trations of the phospholipids in ethanol. 

The aqueous phase used in all the experiments is Milli-Q 
water which has been analysed by light scattering to ensure that 
no extraneous particles are present. The dye used for visualisa- 
tion experiments is a fountain pen ink (Camlin Royal Blue) dis- 
solved in ethanol, which is carefully filtered through a 0.2 fim 
filter paper, to limit the size of particles. After filtering it was 
found that the solid dye particles had a unimodal distribution of 
sizes with mean diameter 220-250 nm. 

2.2. Microfluidic device fabrication procedure 

Here we describe the method to prepare the co-axial flow 
microfluidic device in some detail. The method is easily repro- 
ducible in any laboratory with basic chemicals and instruments. 
An open chamber is prepared with the help of a glass slide 
and strips of double sided tape as side walls. Two Borosilicate 
glass capillary tubes (30 mm and 90 mm in length, 1 mm outer 
diameter (OD) and 0.5 mm inner diameter (ID): GD-1, Nar- 
ishige, Japan) are joined together using using Araldite® (which 
is soluble in chloroform) into a Y-shaped geometry as shown 
in Figure [T] and is inserted through a hole in the tape inside 
the chamber. Liquid Polymer poly dimethyl siloxane (PDMS) 
(Dow corning Sylgard 1 84) is then poured in the chamber and 
left for curing at 90° C for two hours. The glass slides are re- 
moved and the now solid PDMS block is left in chloroform for 
two hours to allow it to swell (degree of swelling is >150%). 
The side glass capillary is now detachable, leaving a channel 
in its place. The main capillary is retracted partly up to the 
Y-j unction, and left inside the PDMS block to serve as a hy- 
drophilic side wall for the outlet stream. A pulled capillary 
(pulled to an ID of 40 /mi, with a glass micropipette puller: 
Narishige PC- 10, Japan), is then inserted through the hollow 
PDMS channel to go past the Y-j unction through the main cap- 
illary. The length of the pulled capillary is such that the en- 
tire tube is inside the PDMS channel. Two blunt needles of 
appropriate size are then inserted through inlet channels up to 
the glass capillary. The PDMS block is then left at room tem- 
perature for the chloroform to evaporate and regain its original 
shape and size, now tightly holding the capillary and the nee- 
dles. Since the flow usually has a small Reynolds number, the 
flow through the side channel attains the fully developed flow in 
a short distance, and can be used to hydrodynamically focus the 



The phospholipids used to synthesize the liposomes in the 
present method are: 1 ,2-dimyristoyl-sn-glycero-3-phosphocholine 
(14:0-DMPC); l,2-dipalmitoyl-sn-glycero-3-phosphocholine (16:0- core flow - Teflon tubin S is used t0 connect the microchannel to 
DPPC) (Avanti Polar Lipids), l,2-dioleoyl-sn-glycero-3-phosphocho t to rin g e 4 pha P al et al.[ |20lT). Two Syringe pumps (NE1000, 
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Figure 1: Top: Schematic diagram of core-annular microfluidic device fabrication sequence. Middle: Final assembled configuration of the entire setup. Bottom: 
Cross section of axial flow focusing and a sample photograph showing the spreading of the dye from the core stream into the annular region. The liposomes are 
expected to be formed in the mixed region. 



New Era pumps Inc.) are used for pumping the fluids inside mi- 
crochannel. 

2.3. Experimental setup for liposome synthesis 

The vesicles synthesis from the solution phase is only pos- 
sible above the phase transition temperature (T g ) of the lipid 
membrane, so it is necessary to maintain a temperature of the 
system above this. When required, the microfluidic device along 
with the long teflon tubings are kept inside an incubator with the 
temperature controlled in the range 27° C to 67° C. 

2.4. Estimation of mean diameter and polydispersity of microflu- 
idic liposomes 

The resultant liposomes are characterized by a single angle 
back-scatter Dynamic Light Scattering (DLS) or Photon Corre- 
lation Spectroscopy (PCS) (Zetasizer nano ZS, Malvern, UK). 
He-Ne laser at 633 nm wavelength and 4 mW power is used and 
the scattered light is detected at an angle of 173°. The size and 
size distribution is calculated using software DTS 6.1 (supplied 
by Malvern). The Polydispersity and the average liposome di- 
ameter are evaluated using cumulant analysis method and the 
weight- weight distribution of diameter is evaluated using multi- 
modal distribution analysis. Three independent experiments are 
carried out for liposome synthesis experiments and for each of 
them, three measurement runs are performed). Ethanol present 
in the dispersion has a large influence on the viscosity of solu- 
tion and a negligible effect on the refractive index; these correc- 
tions are made using the DTS 6.1 software (Malvern) and were 



found to be in good agrement with reported values ( Song and 
Peng) [2008]). 



2.5. Estimation of size and lamellarity of liposomes 

Tunnelling Electron Microscopy (TEM) (JEM-2100F, Jeol) 
is performed to check the lamellarity of the liposomes. An elec- 
tron dense material, 2%w/v phosphotungstic acid (PTA) (SD- 
Fine chemicals India), is used for negative staining . A drop of 
the liposome solution (10 [A) is placed carefully on the TEM 
grid with the help of a micropipette and kept for air drying for 
10 min, followed by placing a drop of negative stain. The sam- 
ple is air dried again before observing at an operating voltage 
of 120 kV. 

Environmental Scanning Electron Microscopy (ESEM) (Quanta- 
250 FEG, FEI) analysis is performed using a Wet-STEM-detector 
and the dispersion is observed at a temperature of 4°C 100% 
humidity, 20 kV accelerating voltage (0.14 nA current), and 
840 Pa pressure. Wet STEM allows observation of liposomes 
keeping its environment fully hydrated i.e. in its native state. 

3. Results 

In the sections below we first report the observations of var- 
ious influences on the size and nature of liposomes formed. The 
co-axial (core-annular) microfluidic device is used for all these 
experiments, unless otherwise stated. Lipid solution in ethanol 
is used as the core fluid stream and water is used in the outer 
stream. 



3 





180 




160 




140 


E 






120 


i- 




0) 



100 


E 




re 


80 


T3 


c 




ns 


60 











40 




20 








A Const Q-annular 
O Const Q-core 
OQ-total varying 







° $ <j> <J> $ (j) 4 



10 20 30 40 50 60 70 80 90 100 110 
Flow rate ratio (4>) 

Figure 2: Effect of the flow rate ratio (outer to inner) and total flow rate on 
the size of liposome formed. The mean diameter is obtained by DLS. Trian- 
gles: <2w (aqueous in the annular region) is kept constant at 100 //1/min and Qe 
(lipidic ethanol in the core region) increases going from right to left. Circles: 
Qe is constant at 10 //1/min and <2w increases from 100 to 1000 /zl/min (going 
from left to right). Diamonds: = 20 is kept constant but the total flow rate Q 
is varied from 21 to 945 //1/min, by increasing both the flow rates. 



3.1. Effect of flow on the mean liposome size 

The ratio of the volumetric flow rates of the aqueous (an- 
nular) and lipids dissolved in ethanol (core) has been shown to 
have a significant effect on the size of liposomes in the rectan- 
gular cross section geometry ( |Jahn et aL) |2007| ). We confirm 
that a similar behaviour is reproducible even in the axial flow 
focussing geometry. In a later section we also provide a plausi- 
ble explanation for this behaviour. 

The flow rate ratio is denoted here as 0, and defined as 







Q, 



annular 



Gc 



Qw 

Qe 



(1) 



where Qannuiar (or gw ) is the outer (water) volumetric flow rate 
and <2 C ore (or Qe) is the inner (ethanol) volumetric flow rate. It 
is seen from Figure [2] that at lower 0, the size decreases with 
increase in 0, and plateaus to a nearly constant value at large 
0. In this study has been varied by two ways, one by keeping 
the inner flow rate constant and other by keeping the outer flow 
rate constant. The liposome sizes have been determined by PCS 
(single angle DLS). 

It is observed that in both the cases (constant gw and con- 
stant Qe), the qualitative behaviour is the same-a decrease fol- 
lowed by a plateau. The large plateau for the liposome sizes 
has not been reported in earlier works of microfluidic synthe- 
sis ( |Jahn et aL) 2007 1 |2010| ), probably owing to the limitation 
of a channel geometry where obtaining a stable flow is difficult 
at high flow rate ratios. Though not shown here, in the axial 
flow focussing we can obtain flow rate ratios as high as 0= 500 
(Qi=5 yul/min and <2w=25()() yul/min). 

Whereas the large plateau of the liposome size is nearly 
equal in both the cases, there is a variation in the size of lipo- 
somes in the two cases for intermediate around 20 < < 30. 
This suggests that there could be other parameters influencing 



the size of liposomes, rather than just the flow-rate ratio, as 
thought previously (Jahn et al. 2007| ). 

The presence of other influences on the size is more clearly 
seen when the flow rate ratio is kept constant, while increas- 
ing the overall flow rate Q = Qw + 2e- We show this in Fig- 
ure [3] where the size of liposomes, at a constant flow rate ratio 
= 20, is plotted against the Reynolds number Re = dUp/ji, 
where U is the superficial velocity U = 4Q/nd 2 , d is the inner 
diameter of the outer channel, p is the mean mass density and [i 
is the viscosity of the fluid (here taken as that of water). At low 
velocities, the size is nearly constant (at a value around 100 nm, 
comparable to the low sizes shown in Figure [2]), whereas at 
higher velocities it approaches a value close to 70 nm, close to 
the large plateau of Figure [2] This variation is also shown 
in Figure [2] (as diamonds at constant = 20). Again this ob- 
servation has not been reported earlier, possibly because in the 
regimes explored in Jah n et al] ( |2007| ) the overall flow rate was 
incidentally varied only in the region where it had no influence. 

The flow rate ratio also specifies the ternary component 
composition (water: ethanol: lipid) in the final collected sample. 
A lower implies a higher fraction of lipid and ethanol in the 
ultimate mixture. The size of the liposomes cannot, therefore, 
be characterised in a simple ternary component diagram pshii 
et al.| 1995) ) which are commonly used for equilibrium systems. 

Increasing the overall flow rate amounts to two obvious in- 
fluences on the local transport phenomena: increase in the local 
shear rate and an increase in the Peclet number (that measures 
the relative importance of axial convection to the radial diffu- 
sion of species) defined in this case as: 



Pe 



Ud 2 



(2) 



where, D is the diffusivity of a solute and L is a characteris- 
tic length along the axial direction. The solute in this case is 
ethanol, which dissolves and diffuses from the core region into 
water in the annular region. 

The influence of shear rate on the dynamics of self assembly 
of lipids into disc-like micelles, and further aggregation of these 
micelles towards larger domain growth and closure to form lipid 
vesicles, can be one possible influence. However, to the best 
of our knowledge this has not been studied, and therefore, we 
summarily leave out exploring this possibility. 

Figure [3] indicates that increase in the Peclet number Pe 
(only U is varied in Eq. ([2])) leads to a smaller sized vesicle. The 
value of Pe for this study is varied from (9(10) to (9(100). A 
large Pe » 1 implies that the two streams hardly mix by diffu- 
sion inside the channel of the device. Mixing of the two streams 
is essential for the formation liposomes. It is only when water 
and ethanol inter-diffuse and the lipids can no longer remain in 
the dispersed form (due to the energetic costs in the presence 
surrounding water molecules) that they self-assemble to form 
vesicles. Figure [3] indicates that Pe values are very large for 
the entire range of flow rates studied, and there is a variation in 
the liposome sizes. This is counter intuitive. Pe » 1 implies 
the flow inside the channel does not lead to (diffusive) mixing 
of the streams and therefore does not induce any self assembly. 
The mixing of the two streams takes place outside the channel 
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Figure 3: Effect of the total flow rate on the size of liposomes. The Reynolds 
number Re is shown as one measure of the total flow rate. Also shown is 
an equivalent axis in terms of the Peclet number Pe. The experiments are 
conducted at a constant = 20, and varying total flow rate Q such that 
Q E = 5,10,...,45,50//l/min. 



as the outlet flow is collected in a cuvette, in which case the 
liposome size should have been uniform for all Pe » 1 . This 
suggests that there could be other phenomena occurring in the 
channel at Pe « 200 that is inducing an increase in the liposome 
size. 

3.2. Influence of the type of Mixing 

Though several variants of ethanol injection methods have 



been employe d ([Batzri and Korn||1973[ |Kremer et al.| |1977[ 
| Wagner et al.| |2002[ |Ishii et al. 1995 ), there is no reported 
analysis of what influences the size of liposomes. Here we 
employed various types of mixing the fluids and observed the 
dependence of the liposome size. Figure [4] shows that the lipo- 
some size is a function of lipid concentration for various ways 
of mixing the lipid/ethanol with the aqueous phase: (a) mixing 
of lipid-ethanol by injection in the aqueous phase under contin- 
uous stirring, maintaining the ratio of volumes at 1:10 (equiv- 
alent to the flow situation of = 10) (b) bulk vigorous mixing 
of lipid-ethanol phase stream with an aqueous phase stream, 
each of them pumped using calibrated syringe pumps. Here the 
flow rate of lipid solution was 10 yul/min and aqueous phase was 
100 yul/min, and (c) Mixing by microfluidic device, here again 
the flow rates were similar as in (b). In each case, the ratios 
of lipid: ethanol: water are kept constant. It is clear from Fig- 
ure [4] that direct mixing by turbulence yields smaller vesicles 
(SUVs), mixing in streams gives slightly larger vesicles and mi- 
crofluidic mixing gives even larger vesicles. The first two meth- 
ods showed a dependence of size on the lipid concentration in 
ethanol but in the case of microfluidic mixing, liposome size 
seems independent of lipid concentration, for even as high as 
80 mg/ml of DMPC. Methods (a) and (b) also show poly dis- 
persed population (not shown here). 

Though the implication of the effect of concentration is not 
evident, we can conclude on the effect of nature of mixing on 
the size of liposomes. The stronger the nature of mixing the 
smaller the size. In order to form larger liposomes, the bilayer 




Figure 5: Sample TEM micrographs of DMPC liposomes showing unilamel- 
lar membrane liposomes. The synthesis conditions are: 10 mg/ml DMPC, 
10 //1/min ethanol flow rate and 200 //1/min water flow rate. 



membrane discs has to grow to large sizes before closing upon 
themselves. However, this growth seems to be hindered by the 
stronger flow, leading to earlier closure and therefore smaller 
liposomes. 

3.3. Liposome Characterisation 

The liposomes formed by the co-axial flow focussing de- 
vice have also been characterised by the usual techniques. The 
particle size distribution by DLS of liposomes is nearly mono- 
disperse with the polydispersity index is around 0.2+0.5. Fig- 
ure|5]shows a Transmission Electron Microscope (TEM) micro- 
graph of DMPC liposomes, clearly displaying the unilamellar 
bilayer membranes using the negative staining method used for 
visualizing. Combined with the near mono-disperse nature of 
the suspension, this method is well suited to carry calculated 
amounts drugs. Figure [6] shows the intact and spherical shaped 
liposome, fairly uniform size distribution of liposome prepared 
with microfluidic synthesis. Figure [6] also shows a comparison 
of particle size distribution obtained with DLS along with that 
obtained by visual analysis of ESEM results, indicating that the 
size distribution obtained in DLS is a reliable measure of the 
actual distribution. 

3. 4. Effect of lipid concentration 

In coacervation techniques, it is well known that liposome 
size and polydispersity increases with lipid concentration ( |Kre-| 
mer et al.) |1977[ | Wagner et al.| |2006| ). We report the effect of 
three different DMPC stock concentration- 10, 20 and 40 mg/ml 
of ethanol using the axial flow-focussing device (Figure [7}. For 
10 and 20 mg/ml of lipids, the concentration does not seem to 
affect the size of liposomes. The concentration of lipids seems 
to have an effect only at large values of 0. At smaller values 
of the concentration of lipids does not significantly influence 
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Figure 4: Three types of mixing methods employed to study the effects of bulk mixing: (a) Bulk Injection — by injecting lipid-ethanol in the aqueous phase under 
continuous stirring; (b) Stream mixing (Syringe Pump) — mixing of lipid-ethanol phase stream with an aqueous phase stream each driven the help of syringe pumps 
and fine Teflon tubing. The flow rate of lipid solution was 10 jA per min and aqueous phase was kept at 100 jA per min. (c) Laminar Microfluidic mixing. The flow 
rates were same as in stream mixing. In all the above methods, the volumetric ratio of ethanol to aqueous phase is 1:10. The plot on the right hand side shows the 
liposome size as a function of concentration of lipids and the type of mixing. 




HV curr I det I temp pressure 1 urn 

20.00 kV 0.14 nA STEM 4.0 °C 840 Pa Quanta FEG 



Figure 6: Environmental Scanning Electron Microscopy (ESEM) of HSPC liposomes: <2e=10 jA/min, <2w=400 fA/mm, 0=40. Right: Comparison of number 
distribution of liposome size obtained by ESEM and DLS. 



the size of the liposomes. This also happens to be the region of 
(p where we see a variation in the size upon changes in 0. This 
again indicates that there are other phenomenon that influences 
self assembly in this region of flow. 

3.5. Effect of membrane composition 

The mean diameter of vesicles is also a function of phos- 
pholipid chemistry and membrane composition as seen from 
Figure [8] To study the effect of cholesterol on liposome size, 
DMPC: Cholesterol in a ratio 9:1 %w/w is dissolved in ethanol, 
so that the final lipid concentration in the stock solution is 10 mg/ml 
We see that addition of cholesterol increases the liposome size. 
Addition of cholesterol increases the bending rigidity modulus 



pared to more flexib l e disk (plain DMPC K b = 0.56 x 10" 19 J 
( [Evans and Rawicz| [T9901 |Zook and Vreeland| |20T(j] )). The 
rigid bilayer disk curves into vesicles by thermal fluctuations 
only when they have grown in size above a critical dimension 
( [Leng et aT] |2003| ). Increased stiffness implies the discs need 
to grow larger (r > r c ) before thermal fluctuation can induce a 
closure. 

Secondly, we alter the lipid chain length and saturation by 
considering other lipids in place of DMPC as shown in Table [T] 
DOPC produces vesicles larger in size than DMPC lipid s. The 
bending modulus for DOPC is K b = 0.85 x 10" 19 J ( |Rawicz| 



et al. 2000), which is larger than that for pure DMPC. Similarly 
DPPC which has a smaller bending modulus K b = 0.4 x 10" 19 J 



Kb and area compressibility modulus K a ( [Tierney et al.| [2005 ). 
Increase in K b increases stiffness (DM PC+Cholesterol Kb = 
4.2 x 10" 19 J ([Evans and Rawicz[ [1990])) thus it becomes diffi- 



cult for a stiffer bilayer disk to bend and form vesicle as com- 



(Lindahl and Edholm, 2000) results in smaller sized liposomes. 
In general, it is expected that the average diameter of vesicles 
increases with increase the bending rigidity (Helfri ch]|1973| ). 
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Table 1 : Chain length and saturated bonds for various lipids employed in this work. 
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Figure 7: Effect of phospholipid concentration (in ethanol) on the liposome 
size. At lower 0, size is insensitive to the concentration whereas at higher 0, 
the higher concentration of 40 mg/ml produces larger sized liposomes. 



Figure 8: Liposome size as a function of flow rate ratio and membrane com- 
position. 



4. Discussion 

4.1. Hypothesis for Liposome Size Variation 

The results of the previous section can be summarised as 
follows (i) there are two regimes of liposome formation: one 
in which the size depends on the flow parameters but does not 
depend on the concentration of lipids, and other in which the 
size is not dependent on flow parameters but is influenced by 
concentration of lipids, (ii) The dependence on the lipid chem- 
istry and membrane composition is in line with what is expected 
from a general understanding of the bilayer growth and closure 
to for vesicles. Therefore, it is the former observation that needs 
an explanation. 

We hypothesise that the two regimes correspond to (a) Con- 
vective instability mixing: Instability induced convective mix- 
ing within the channel, and (b) Droplet mixing: Mixing of the 
streams outside the channel when drops of the outlet are col- 
lected in a cuvette. Core-annular flows of two miscible fluids 
where there is a viscosity stratification leads to an instability at 
large Peclet numbers (Pe » 1) and small Reynolds numbers 
(Re = dU pip). This was numerically shown by Ranganathan 



and Govindarajan (2001 ) for channel flows and by |Selvam et al.| 
(2007]) for core-annular flows. In the present case, there is a 



viscosity variation in a non-monotonic manner, unlike that con- 
sidered in Selvam et al. ( 2007 ) where the viscosity variation is 
monotonic. |d'Qlce et al.| (|2008) have observed certain insta- 
bilities experimentally which show pearl and mushroom-like 
appearance when there is a monotonic increase in the viscos- 
ity from the core region to the annular region. To the best of 
our knowledge such instabilities have not been reported in the 
water-ethanol system. Nevertheless, we may expect a similar 
behaviour. 

The reason for expecting the droplet mixing regime is that 
when the Pe number is large and when there is no convective 
instability, the streams will only mix at the outlet. This is sup- 
ported by additional facts that the dependence on the concen- 
tration of lipids is similar to that seen when two streams are 
mixed in the various types of ethanol injection methods shown 
in Figure]?] 

At present we do not have a theory to show the possible 
existence of either of the above claims for the present system. 
However, it is possible to test the consequences by carrying out 
visualisation experiments on the same system. In the follow- 
ing we report results of experiments carried out in a similar 
system — the lipids in ethanol being replaced by a dye — which 
us allows to visualise the nature of flow inside the device. 
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Figure 9: Viscosity (Song and Peng 2008) and density (Perr y et al.] |1997) 
variation of ethanol- water system at different temperatures. 



Figure 10: Validation of the width determination. Here the flow rate ratio is 
kept constant at 0=50, increasing the total flow rate. 



4.2. Flow Instability due to Viscosity Stratification 

Core-annular instability is known to occur in immiscible 
flows. It was shown in Selvam et al.| p007) a similar instability 
also occurs in miscible fluids, and presence of a small amounts 
of diffusivity (Pe » 1) enhances this instability, particularly at 
small Reynolds Numbers. Though the model is not identical to 
the present system, we draw similarities in the nature of flow 
(core-annular) and the presence of viscosity stratification. In 
the present case 2 < Re < 42 and 16 < Pe < 507, based on 



water-ethanol diffusivity of D = 1.28 x 10 9 (Zhang et al. 

[2006] ). 



The variation of viscosity of water-ethanol system ( |Song| 
and Peng] |2008| ) is shown in Figure [9] The maximum viscos- 



ity of water-ethanol mixture at 20°C is 2.68 mPa.S at 0.7 mole 
fraction of water; which is around three times the viscosity of 
water. |d'Qlce et al.| (|2008|) observed pearl and mushroom like 



patterns of the instability in the case of a water-natrosol system 
where the core fluid was of a lower viscosity, and the viscosity 
increases monotonically to about 25 times the core fluid viscos- 
ity. |Selvam et al.| ( |2007"] ) also reported instabilities to be present 
in the opposite case of viscosity decreasing from the core to the 
annular region. In the present case, we have both the situations. 
The viscosity is maximum in a small mixed region in between 
(at about 40%w/w of ethanol), and decreases to about the same 
value (1 cp) towards core and the periphery. 

To study the nature of flow, we carry out the experiments in 
the absence of lipids, thereby not including the additional com- 
plexity of liposome formation. Instead, we introduce a dye — 
filtered fountain pen ink — to enable visualisation. Experiments 
were conducted at various flow rate ratios identical to that em- 
ployed in Figure [2] Since the dye is introduced in the core fluid, 
we expect that in absence of diffusion of dye, or any other con- 
vective instability, the width w of the core (coloured) region in 
a fully developed flow would be 



(3) 



This can be obtained by integrating the fully developed Poiseuille 
flow parabolic velocity profile (assuming uniform viscosity) for 
the core and annular-fluid volumetric flow rates, and in the limit 
of w <?c d. We validate the experimental method of determining 
the width of the dye-region (core fluid) which results in the data 
shown in Figure [TO] We choose a large value of = 50, and 
vary the total flow rate. Here, 17 < Pe < 507 and 2 < Re < 64. 
As expected, the diffusive mixing should be minimal for large 
Pe, and the observed width is in agreement with the expected 
value from Eq. ([3]). 

We have considered a similar set of flow rates as used in the 
liposome synthesis in Figure [2] In Figure 1 1 we have consid- 
ered three variations: (a) Constant flow rate ratio cp = 20 and 
cp = 50, (b) Varying cp keeping outer flow rate constant, and (c) 
Varying (p keeping the inner flow rate constant. At a large flow 
rate ratio of (p = 50, we see that the width of the core region 
remains the same. 



In Figure[TT](c) we see an unexpected behaviour. The width 
of the dye region increases with increase in the inner flow rate. 
We expect that this will decrease (p from about 100 to about 10. 
Though this implies an increase in w from Eq. ([3}, the increase 
is much more than what is theoretically expected, as shown 



in Figure 12 In this set of experiments the Peclet number is 
high at Pe « 35, and increases only marginally from Pe=33 to 
Pe=36, so this increase in the width cannot be explained even 
by Taylor diffusion of ethanol into water. Moreover, the dye 
particles being larger than ethanol, the Peclet number for dye 
particles is even larger at Ped ye > 10 3 . Therefore, the explana- 
tion of the increase in the observed width cannot be given based 
on mere mass conservation or diffusion, and we conclude that 
a gentle convective mixing of the fluids takes place due to the 
growth of an unstable mode. This is a gentle mixing because by 
nature of the flow, once the fluids mix at low Re the viscosity 
becomes uniform, which stabilises the flow. 

A similar situation is seen when the outer flow rate is varied 



keeping the inner flow rate constant, as seen in Figure 1 1 (d) 
Here too the Peclet number is high, 36 < Pe < 335. While 
we can conclude the existence of a convective instability that 
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leads to mixing of the fluids, it is difficult to correlate this to the 
theoretically expected instability value of the Pe or the Re, as 



the specific model considered in |Selvamet a l. ( 2007 ) is different 
from the present system. 

The mixing of the fluids can be schematically explained 
with the help of Figure [14] In the stable case with high Pe, 



the fluid packets move in a laminar fashion, and the diffusion of 
the particles out of the packet is negligible. Whereas, when the 
viscosity stratification leads to an instability, the flow does not 
remain laminar — it mixes the packets of fluids. Here too since 
the Pe is large, the diffusion of particles out of the packet is 
negligible on the length scale of the channel diameter, however 
when seen from a macroscopic perspective, the dye particles 
occupy a wider region in space as the packets have themselves 
been dispersed by mixing. To contrast this, a diffusive mixing 
occurs at a low Peclet and laminar flow, where packets of fluid 
exchange particles from other packets, purely by diffusion. 

4.3. Droplet mixing regime 

From the above observations we conclude that not all high 
Pe and low Re cases leads to an instability. There are regions 
where there is an agreement of the theoretical width with what 



is observed. This is summarised in Figure 13 where we see 
that for very small w there is close agreement, whereas it begins 
to deviate at w/d around 0.3. In the regime of agreement, the 
width of the core flow remains nearly constant throughout the 
length of the channel. At the outlet of the channel, the flow 
form droplets (the overall flow rates being small). Since the two 
streams remain unmixed inside the channel, even while forming 
the drops we observe two visually distinct regions. The fluids 
only mix when the droplets are gathered in the outlet reservoir. 

In Figure [l3j the correlation between the liposome size and 
the width of the dye can be observed. Though these are two dif- 
ferent experiments, they have been carried out under the same 
flow rate conditions. The liposome size is plotted against the 
observed width, each obtained under the same flow rate val- 
ues. This plot can be used to infer the nature of flow that would 
have occurred in the liposome synthesis experiment. When the 
agreement of the observed width with the theoretical width is 
good the liposome size is on the lower side, implying the li- 
posome size is governed more by the droplet mixing regime. 
Moreover, when the observed width deviates from the theory, 
the liposome size is the largest. This implies that it is the inter- 
nal convective mixing that is leading to larger size. 

4.4. Hypothesis for liposome size variation 

Combining the observations of the above experiments with 
the size of liposomes, we conclude that the liposomes are smaller 
when the streams mix in the droplet regime, and are larger and 
when the streams mix by convection inside the channel. The 
size of liposomes formed in the droplet mixing regime are also 
comparable with those obtained in the ethanol injection and 
stream mixing regimes shown in Figure [4] 

We hypothesise that the size of the liposomes depends on 
the nature of convective mixing. More gentle and less turbu- 
lent the mixing, the larger the size. In the cases of droplet, 
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Figure 1 1 : Visualisation of mixing of a blue dyed ethanol (without the lipids) 
with water, in four typical cases of flow rates employed for liposome synthesis 
in Figures [2] and [3] (a) Low = 20 kept constant, and total flow rate Q is 
varied such that Qe = 5, 10, . . . ,45, 50 //1/min, going from left to right. The 
photographs shows wider region of mixing at lower Q. (b) High = 50 kept 
constant, and total flow rate Q is varied such that Qe = 2, 4, . . . , 16, &30 //1/min, 
going from left to right. The photographs show no mixing of the dye. (c) Con- 
stant <2w = 100 //1/min and increasing Qe = 1,2, ...,8,9 //1/min showing 
strong mixing for low (right most), (d) Constant Qe = 5 //1/min and increas- 
ing <2w = 100,200, . . . , 1000 //1/min. showing strong mixing at low (left 
most). 



ethanol injection, and stream mixing, the flow is locally vigor- 
ous, whereas in the viscosity stratification led instability mixes 
the fluids rather gently. 

As water mixes with ethanol, the lipids self-assemble when 
it is no longer favorable to remain in a dispersed form. They be- 



gin to assemble as disc-like bilayer membrane micelles ( Leng 
et al. , 2003) which grow in size. When the edge energy at the 
periphery of the discs become comparable to thermal fluctua- 
tions induced bending energy of the membrane, the discs close 
upon to form vesicles ( Helfrich| |1973| ). Left to itself all the vesi- 
cles formed would be about the same size as determined by this 
balance. However, convective motion in the length scales of the 
discs can induce closure before the membrane grows to larger 
possible dimensions. It is known that the smallest vesicles are 
induced by the action of ultrasonic waves, that essentially intro- 
duce convective motion in small length scales by spontaneous 
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Figure 12: Comparison of the observed spread of the dyed ethanol with near w/d (Observed) 

the outlet of the channel, along with the theoretical approximation of the width 
obtained from Eq. J3j. 



vapourisation into bubbles ( [Richardson et al.||2007| |Yamaguchi 
[eTaTl[2009] ). 

The presence of ethanol is to essentially reduce the interfa- 
cial tension along the edges. A reduced tension will lead to a 
larger sized liposome, since the thermally induced bending of 
the membrane takes place at a larger size of the disc. 

5. Conclusion 



Agreement line 

Theoretical Width 
A Liposome Size 

Figure 13: Correlation of the mixing observed visually in Figure [12] with the 
liposome size obtained under identical flow conditions. Liposome size are 
smaller when the theoretical width matches with the observed (droplet mix- 
ing). Liposomes of largest size are formed when there is a large disagreement 
in the widths (convective mixing). 



Two distinct but connected aspects of microfluidic synthesis 
of liposomes were shown. The general method employed for li- 
posome synthesis is a variant of the ethanol injection technique 
in which lipids dissolved in ethanol is contacted with water, 
leading to the assembly of lipids to closed vesicular structures 
in the presence of excess water. 

The first aspect is a detailed demonstration of a simple de- 
vice capable to synthesise unilamellar liposomes of nearly uni- 
form size in a co-axial core-annular flow geometry. Detailed 
observations of the effect of various flow parameters, lipid chem- 
istry and concentration were made. The chief conclusion from 
these observations are that the size and formation of liposomes 
cannot be characterised by ternary component diagrams, or just 
the flow rate ratio as done in the past; and that the size strongly 
depends on the nature flow. 

The second and more significant aspect is on throwing some 
light on the hitherto unknown mechanism of liposome size se- 
lection in ethanol injection techniques and by implication in a 
microfluidic setup such as the one employed here. We showed 
that the size depends strongly on the nature of micro-convective 
mixing. 

In the case of the microfluidic set up, there are two regimes. 
One in which a convective hydrodynamic instability due to vis- 
cosity stratification leads to a gentle mixing and a larger size 
of liposomes. In the other regime the streams do not mix in- 
side the channel, but mix as droplets formed at the outlet when 
collected in a reservoir. Vigorous mixing leads to smaller lipo- 
somes. Though we do not present a detailed theory for the exact 



value of size selection, it is argued that it is the length scales of 
flow (mixing) that determine the size of liposomes. 

Further studies into the detailed structure of the flow in mi- 
crofluidic and other ethanol-injection methods can establish a 
more accurate connection between the liposome size and the 
length scales of flow variation. The discovery of this impor- 
tant effect should help in formulating models that will be able 
to predict liposome size and lamellarity in many coacervation 
methods. 
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